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C " 3 Low angular momentum accretion flows can have standing and oscillating shock waves. 

, We study the region of the parameter space in which multiple sonic points occur in vis- 

C " D ■ cous flows in presence of various cooling effects such as bremsstrahlung and Comptoniza- 

tion. We also quantify the parameter space in which shocks are steady or oscillating. We 
find that cooling induces effects opposite to heating by viscosity even in modifying the 
topology of the solutions, though one can never be exactly balanced by the other due to 
^ ■ their dissimilar dependence on dynamic and thermodynamic parameters. We show that 

' beyond a critical value of cooling, the flow ceases to contain a shock wave. 
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1. Introduction 

Cooling and heating processes play an important role in studying the accretion 
disks around compact objects. In this paper, we consider both viscous lieating and 
bremsstrahlung cooling as energy dissipative processes. In Chakrabarti and Das^, 
the problem of accretion and winds with small angular momentum was solved in 
a very comprehensive way when viscous heating is included. The entire parameter 
space was scanned and separated in terms of flows with and without shocks. A 
major scope for further work is to find and separate the parameter space which 
allows shock formations when cooling effects are also included. This will allow us to 
study flows with high accretion rate as well. The question arises because viscosity 
transports angular momentum and increases the possibility of shock formation at 
a larger distance from the black hole. However, cooling reduces the post-shock 
pressure and therefore the possibility of shock formation. Various cooling processes 
will and should change the parameter space in which shocks form. 

In Chakrabarti^, some of the effects of cooling was included. Assuming that 
cooling at a given radius of an accretion disk is proportional to the heating, all pos- 
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sible topologies were presented. While cooling will and should depend on heating, 

there is no reason to believe that this proportionality is independent of radial dis- 
tance. It is therefore necessary to re- visit the problem with explicit form of heating 
and cooling included. 

In the present Paper, we make these important extensions of the previous work 
and show that shocks are still possible in a very large part of the parameter space. 
We find topologies of solutions which are similar to what was found for spiral shock 
study^, but otherwise new for axisymmetric situation. We also find new topologies 
which were not anticipated before. In more recent years, it has become evident that 
the standing shocks may be very important in explaining the spectral properties of 
black hole candidates^ as the post-shock region behaves as the boundary layer where 
accreting matter dissipates its thermal energy and generates hard X-ray by inverse 
Comptonization. This region is also foTind to be responsible to produce relativistic 
outflows. Furthermore, numerical simulations indicated that the shocks may be 
oscillating at nearby regions of the parameter space in presence of cooling effects^ 
and the shock oscillations can also explain intricate properties of quasi-periodic 
oscillations^. Recent observations do support the presence of sub-Keplerian flows 
in accretion disks'""*. 

The present work is done around a Schwarzschild black hole by using pseudo- 
Newtonian potential®. We use a similar viscosity prescription as in Chakrabarti^. 
In §2, we present model equations which included both heating and cooling effects. 
In §3, we perform the sonic point analysis. In §5, we study the nature of the sonic 
points and how they vary with flow parameters. In §6, we study the global solution 
topologies with heating and cooling effects. In §7, we classify the parameter space 
in terms of whether shocks will form or not and how it depends on flow parameters. 
In §8, we briefly report how even the super-Keplerian flows may also be transonic. 
Finally, in §9, we discuss the importance of these solutions and make concluding 
remarks. 

2. Model Equations When Cooling Effects are Included 

As far as the cooling processes are concerned, they could be due to various physical 
reasons, such as the thermal and the non-thermal bremsstrahlung, synchrotron, 
Comptonization etc. For simplicity, we assume that the Comptonization enhances 
the injected photon intensity due to bremsstrahlung by a factor of C, which can take 
any value from 1 to ~ few xlOO depending on the availability of soft photons^. In 
other words, we use C as a parameter to represent the net cooling. 

We start with a thin, axisymmetric, rotating viscous accretion flow around a 
Schwarzschild black hole. The space-time geometry around a non-rotating black hole 
can be satisfactorily described by the pseudo-Newtonian potential® and is given by 
g{x) = — 2(a-i) ' where, x is the radial distance in dimensionless unit. In the steady 
state, the dimensionless hydrodynamic equations that govern the infalling matter 
are the foUowings^. 



2, 2008 13:30 WSPC/INSTRUCTION FILE ws-ijmpd 



Properties of Accretion Shocks in Viscous Flows with Cooling Effects 3 



(a) Radial momentum equation : 
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apart from the geometric constant. 

(c) Angular momentum conservation equation : 
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(d) The entropy generation equation : 




(x^W^^) = 
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s,?r^ = Q+-Q 



(Id) 



Here, •& is the radial velocity and X{x) is the specific angular momentum of the 
flow. The distances, velocities and masses are made dimensionless by using = 
2GMbh /c^ , the Schwarzschild radius, c, the velocity of light and Mbh, the mass 
of the black hole, respectively. Here S and Wxcf, arc the vertically integrated density 
and viscous stress, s is the entropy density of the flow, T is the local temperature, 
(3+ and Q~ are the heat gained and lost by the flow, and M is the mass accretion 
rate. We assume that the accretion flow is in hydrostatic equilibrium in the vertical 
direction and the vertical velocity component is much smaller compared to the radial 
component. With this assumption, the local disk height is obtained by equating the 
pressure gradient force in the vertical direction with the gravitational force. The 
half thickness of the disk is then given by. 



adiabatic index, pressure and density respectively. 

In the present paper, we follow a similar viscosity prescription as given in 
Chakrabarti^ where W^^ = — anH is used. This prescription ensures that viscous 
stress remain continuous at the flow discontinuity (shock) in presence of significant 
radial motion of the accretion flow. 

3. Sonic Point Analysis 

A black hole accretes matter either from its binary companion or from the sur- 
rounding ambient medium. This matter starts with a negligible radial velocity at 
the outer edge of the disk. But it enters into the black hole with a velocity of light^°. 
This ensures that inside the accretion disk there must be at least one point where 



h{x) = ax 



'/'{x-1), (2) 



where, the sound speed is defined as a 



= Vt^TP) where 7, P and p being the 
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the radial velocity exactly matches with the sound speed. This point is known as 

the sonic point. Accretion flow which passes through a shock wave the must cross 
sonic points at least twice. In other words, the flow may be called multi-transonic. 
For the accretion flow, entropy equation (Id) can be simplified as. 



1? 



7-1 



IdP 

p dx 



dp 
p^ dx 



ph 



C-H. 



and then H{= Q~^/ph) takes the form. 



H = Ax{ga^ + 7??^) 



dQ. 

dx 



(3) 



(4) 



where, A = — an-y and g = -^t^- Here, Vl{x) is the angular velocity of the accreting 
matter at the radial distance x, n is the polytropic index (n = ^^j_^^ ), /„ and In+i 
come form the vertically averaged density and pressure. The general expression of 
In is given by^^: 

I - (5) 
^""(2n + l)!- 

Simultaneously with viscous heating, we use Comptonization of the 
bremsstrahlung radiation as the physical cooling process. The following analysis 
is carried out with non-dimensional cooling term C(= Q~ / ph) as 



1?X3/2(X-1)' 



with 



B = 1.4 X 10-2^ (^) 
\ 2k J 



1/2 M 



2-Kml2GcMBH' 



(6a) 



(66) 



p is the mean molecular weight, nip is the mass of the proton and k is the Boltzmann 
constant respectively. 

3.1. Sonic Point Condition 

From equations (la-Id), the sonic point conditions are derived following the general 
procedure^^ and are given by. 



<M 
dx 



N 



where the numerator N is 



N = -- 
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and the denominator D is, 

^_ 2aH (7 + 1)193 



(86) 



(7-1) (7-1) 

At the sonic point, both the numerator and denominator vanishes. For D = 0, 
one can get the expression for the Mach Number M{xc) at the sonic point which 
is given by. 



-nib - V'ml- '^ragmc 
2ma 



7+1 



for an — > 0, 



where, 



= -Aanl^in - 1)(2<7 - 1) - 7(7 + 1), 
mfo = 27 - 2Aa„sf7(7 - l){g - 1), 
rUc = Aa„g'^{'y- 1). 



(9) 

(10a) 
(106) 
(10c) 



In order that the Mach number to be physically acceptable, we do not use the 

negative sign within the square root. In the weak viscosity hmit, Mach number at 
the sonic point reduces to the result as obtained in Chakrabarti^^. 

Setting N = 0, we get an algebraic equation for sound speed at the sonic point 
which is given by. 



where, 
A=- 



F{£c, Ac, Xc) = Aa'^ix) + Ba^{x) + Ca^ix) +V = 0, 



anAjg + jM'^y ^ an^(5a; - 3){ff + 7M^} ^ M'^{5x-3) 
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(11) 

(12a) 

(126) 
(12c) 
(12d) 



We solve the above quadratic equation to obtain the sound speed at the sonic 
point. Das, Chattopadhyay and Chakrabarti^'' found that depending on a given set 
of initial parameters accretion flow may have a maximum of four sonic points where 
one of the sonic points always lies inside the black hole horizon for non-dissipative 
accretion flow"'^'''"'^^. In our present study we also expect the similar result if we are 
interested in the weak viscosity and weak cooling limit. Below, we study the nature 
of the sonic points. 
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4. Nature of the sonic points 

Accreting matter begins subsonically from tlie outer edge of tlie disk and becomes 

supersonic after passing through the sonic point before entering into black hole. 
A flow may contain multiple sonic points depending on the initial set of input 
parameters. Nature of sonic points depends on the value of the velocity gradients 
(d'd/dx) at the sonic points. At each sonic point {M/dx) has distinctly two different 
values. If both the velocity gradients are real and opposite signs, the sonic points is 
saddle type and one is used for accretion and the other used for winds. Nodal type 
sonic point belongs to the class when the derivatives are real and of the same sign. 
When the derivatives arc complex, the sonic point is of spiral type. For a standing 
shock to form, an accretion flow must have more than one saddle type sonic points. 

In Figure la, we plot the variation of specific angular momentum (Ac) as a 
function of the logarithmic sonic point location {xc) for a given viscosity parameter 
(an = 0.1), and dimensionlcss accretion rate (made dimensionless by the Eddington 
rate) to = 1.0 and a given specific energy = 0.013 at the sonic points. The variable 
used in cooling efficiency factor ( (Eqn. 6a). From the bottom curve to the top, 
C = 1, 20, 40, 60, 80 and 100 respectively. In Fig. lb, the same curve is shown 
for an — 0.5, other parameters remaining the same. The long-dashed curve at the 
top represents the Keplerian angular momentum distribution which is completely 
independent of the initial flow parameters and depends only on the geometry. Solid 
part of the curves represents the saddle type sonic points, dotted curve represents 
the nodal type sonic points and the short-dashed curves are for the spiral type 
sonic points. Here, at a higher viscosity, the number of sonic points becomes three 
even with very low angular momentum. For no Comptonization (lowermost curve), 
the viscous heating is so strong that only outermost sonic point (solid part of the 
curve at large radius) exists. Only a large degree of cooling can compensate for the 
viscous heating to bring back the innermost sonic point. 

First notice that the sonic points occur at angular momentum below Keplerian 
value. For lower values of cooling at the sonic point, an accretion flow contains all 
the three types of sonic points in a systematic order: saddle — nodal — spiral — 
nodal — saddle for monotonic increase of location of sonic points. With the increase 
of C the region of spiral type sonic points gradually decreases and finally replaced 
by the nodal type sonic points though multiple sonic points still exist. For further 
increase of C all the nodal type sonic points also disappear and replaced by saddle 
type sonic points. Note that the angular momentum is always sub-Keplcrian. Later, 
we shall show that when the cooing is very strong, sonic points will form even for 
super-Keplerian flows. 

We continue our investigation of the transonic nature of the flow by showing 
in Fig. 2 a series of curves where the specific energy at the sonic point is changed 
(marked). The long dashed curve is the Keplerian distribution as before. Meanings 
of solid, dashed and dotted curves are the same as before. For negative energies 
there are two sonic points, the inner one is saddle type (shown in solid curve) and 
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Fig. l(a-b): Variation of specific angular momentum (Ac) as a function of the logarithmic 

sonic point location (xc) for the viscosity parameter (a) an = 0.1 (left panel) and (b) 
an ~ 0.5 (right panel). Dimensionless accretion rate m = 1.0 and energy at the sonic 
point £c = 0.0013 are chosen. Long-dashed-dotted curve in the upper part is the Keplerian 
angular momentum distribution. Solid curves represent the saddle type sonic points, dotted 
curves represent the nodal type sonic points and the short-dashed curves are for the spiral 
type sonic points. The Comptonization cooling factor is 1 (bottom curve), 20, 40, 60, 
80 and 100 (top curve) respectively. Cleaxly, higher cooling and higher viscosity remove 
the outer sonics. Eventually the disk becomes a Keplerian disk passing through the inner 
sonic point. 



the dashed curve is the spiral type. For each energy, two curves are drawn. The 
thick curve is drawn when both the heating and cooling are included while the thin 
curve is drawn when only the heating is included. The motivation is to impress that 
the character of a solution can be changed when cooling is included. For instance 
the solution for £c = 0.019 with heating and cooling has no spiral or nodal sonic 
points. But when the cooling is turned off, the saddle type point becomes nodal 
type. 

In our further study of the nature of the sonic points we draw in Figs. 3(a-b) 
the variation of energy at the sonic points with a cooling factor. In Fig. 3a, low 
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Fig. 2: Variation of tlie specific angular momentum (Ac) at the sonic point (xc) as a 
function of the specific energy {£c) of the flow. For each energy, two curves are drawn. 
The thick curves axe drawn for the cases when both the heating and cooling are included 
while the thins curves are drawn when only the heating is included. For negative energies 
there are only two sonic points. 



viscosity (an = 0.1) is used, while in Fig. 3b, high viscosity (an = 0.5) is used. 
Other parameters are: A = 1.65 and m ~ I. The curves from bottom to the top are 
for the cooHng factor ( = 1, 20, 40, 60, 80, 100 and 120. Notice that three sonic points 
occur only when the specific energy is positive, i.e., for sufficiently hot Keplerian 
disks or sub-Keplerian flows. For high viscosity outer sonic points almost disappear 
from regions close to the black hole. 

In Fig. 4, we show how the number of sonic points, reduced due to viscous 
heating process, is recovered back with the introduction of cooling. The curves are 
drawn, from the bottom to the top, for C = 1; 20, 60, 80 100 and 140 respectively. 
Solid and dotted curves are for saddle and nodal type sonic points respectively. 
Other parameters are an = 0.1 and to = 1. Polytropic index n = 1.75 and specific 
energy 8 = 0.0018. The long dashed curve gives the boundary of the angular mo- 
mentum Xcrit = 1.1733 and cooling factor ( = 4.2 below which there are no triple 
sonic points, i.e., no shocks in steady flows. 
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log(xJ log{xJ 

Fig. 3(a-b): Variation of the specific energy (£c) at the sonic point (xc) as a function of 
the coohng rate of the flow, (a) an = 0.1 (left panel) and (b) an = 0.5 (right panel). For 
high viscosity, the outer sonic points almost disappear. 

5. Global Solution Topology 

Study of shock properties require multi-transonic accretion flow. Accretion flow 
passes through two different saddle type sonic points and discontinuous jump of the 
flow variable joins these two different branches — one passes through the inner sonic 
point and the other passes through the outer sonic point. This discontinuous jump 
usually known as standing shock transition. In this paper we discuss the nature of 
solution topology in presence of viscous heating and bremsstrahlung cooling. 

In Fig. 5a, we show how the solution topologies change with cooling. Here we 
plotted the Mach number as a function of the logarithmic radial distance for ^ = 1, 
10, 25 and 50 respectively. Other chosen parameters are: Xin = 2.71, Ai„ = 1.68, 
an = 0.01, m = 0.5 respectively. Wc note that the topology opens up to allow flows 
to enter into black holes through the inner sonic points. 

In Fig. 5b, we show six panels in which we assume higher accretion rate and 
higher inner sonic point locations. The parameters are: C = Ij 10, 25, 33.1, 50 and 
70 respectively. Other parameters are: Xin = 3.5, Ai„ = 1.68, an = 0.01, rh = 0.2. 
Here, as the cooling is increased, the topologies open up similarly, but the route to 
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log(xJ 

Fig. 4: Recovery of outer sonic point points as cooling processes are introduced. Solid 

and dotted curves arc for saddle and nodal type sonic points respectively. Parameters are 
marked on the plot. Long dashed curve at Ac = 1.1733 gives a boundary below which 
there are no triple sonic points. 

opening up is difiFerent. For instance, the solution in the fourth panel, with ( = 33.1 
is completely new and intriguing. In this case the flow has the potential to join with 
a Kcplcrian disk far away (with low Mach number), while at the same time, it also 
passes through the inner sonic point. But it has multi- valued solution: there are 
two Mach numbers at a given radial distance in some region. We conjecture that 
this type of solution should be unstable and would cause non-steady accretion. 

Many of these topologies show 'multiple crossings' very similar to what was 
found in the study of spiral shocks^. Actually, since entropy is changing along the 
flow, the two dimensional nature of the plots in Fig. 5 is slightly misleading. In Fig. 
6, we show a 3-dimensional plot of the second panel {( = 10) of Fig. 5a in which 
the specific entropy is also plotted as it varies along the flow. In particular, this 
diagram shows that at the true sonic point, the specific entropy is identical in both 
the incoming (solid) and outgoing (dotted) branches. But at the 'intersection' (at 
around M ~ 0.5) marked by two squares, the entropies are completely different. 
Thus, there is no 'sonic point' around M ~ 0.5. 
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Fig. 5: Solution topologies in presence of heating and cooling. In (a), the parameters are 
^ = 1, 10, 25 and 50 respectively. Other parameters arc Xj„ = 2.71, A = 1.68, a = 0.01, 
rh = 0.5. In (b), the parameters are C = 1: lOi 25, 33.1, 50 and 70 respectively. Other 
parameters are; Xin = 3.5, A = 1.68, a = 0.01, rh = 2.0. 

6. Parameter Space Description 

It is useful to study the global behaviour of the accretion solutions. For this we 
classify the parameter space spanned by the specific angular momentum and energy 
at the sonic point. Figure 7 shows the parameter spaces in which the solution 
passing through the inner sonic point contains a closed spiraling topology as in 
panel 1 of Fig. 5a. This means that whether the Rankine-Hugoniot relation is 
satisfied or not, a shock could form in this region. The shock will be stationary 
if the Rankine-Hugoniot relation is satisfied^^ and will be oscillating if the relation 
is not satisfied. The solid, dashed, dot-dashed and dot-long-dashed regions are for 
cooling parameter 1, 10, 30, and 50 respectively. As the cooling is increased, the 
region shrinks and becomes smaller and smaller. This indicates that there exists 
a critical cooling parameter, beyond which a flow will cease to have three sonic 
points. Other parameters used are an = 0.01 and rh = 0.1. When the viscosity and 
cooling are reduced to zero, this region merges exactly to the corresponding region 
in C89. 
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Entropy 




Fig. 6: Three-dimensional plot of the second panel (C = 10) of Fig. 5a in which the 

specific entropy is also plotted as it varies along the flow. The boxes represent the pseudo- 
intersection point of that panel at around M ^ 0.5. The two branches have different 
entropies. 

In Fig. 8(a-c), we present a few complete solutions which are drawn with the 
parameters at Ac = 1.7, an = 0.05, m = 0.2 and ( = 5 and only the inner 
sonic point is varied: (a) Xin = 2.545, (b) Xin = 2.55 and (c) Xin = 2.555. The 
corresponding shock locations are (a) Xg = 48.199, (b) Xg = 27.8854 and (c) Xg = 
18.6445 respectively. Vertical dashed lines show the shock transitions which connect 
two solutions, one passing through the inner sonic point and the other passing 
through the outer sonic point. The spiral loop through the inner sonic point rapidly 
shrinks with the increase in the sonic point location. The shock location also comes 
closer. This shows that even though the shock location may change by orders of 
magnitude, the inner sonic point virtually remains at the same place. 

It is instructive to know the sub-division of the parameter space in terms of the 
topologies of the solutions. All the topologies seen in Paper I are also present in 
this case, but a new topologies occur (see. Fig. 5b) when the cooling is especially 
strong. Figure 9 shows the sub-division of the parameter space which arc marked 
and the corresponding topologies are shown in the bottom left. When the cooling is 
very strong, the curve ABC shows further sub-divisions and a new topology shown 
in the box C13 occur. The number of loops in the inflow may increase depending 
on the cooling. The dotted curves indicate that the Figure is drawn for different 
(high) cooling factor (■ The regions marked S, OS, Cll, C12, OAC, /*, O* 
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Fig. 7: Region of the parameter space in which the solution passing through the inner sonic 
point contains a closed spiraJing topology. The solid, dashed, dot-dashed and dot-long- 
dashed regions axe for cooling parameter 1, 10, 30, and 50 respectively. Other parameters 
are an = 0.01 and rh = 0.1. 

produce topologies which produce standing shocks in accretion, oscillating shocks 
in accretion, region which produces one type of closed topology (clockwise turn), 
region which produces the other type of closed topology (anti-clockwise turn) , region 
which produces one open and the other closed topology, region which produces open 
solutions passing only through the outer saddle type sonic point and the region 
which produces open solutions passing only through the inner saddle type sonic 
point respectively. The region producing the new dotted topology C13 is very close 
to the curve ABC and can be discernible only when cooling is strong. 

It would be of interest, to concentrate on the modification of the parameter space 
for shock formation when cooling is enhanced. In Figs. 10 (a-b) this is shown. In Fig. 
10a, an = 0.01 is chosen and the cooling parameters are ( = 0.01, (dot-dashed) 0.1 
(long-dashed) and 1 (solid). We note that the region of the parameter space shifts 
to include negative energy regions as well (For instance, for au = and C = 0) the 
parameter space contains only positive energy.). In Fig. 10b, viscosity parameter is 
increased to an = 0.05. This causes a shrinkage in the parameter space. The cooling 
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log(x) 

Fig. 8(a-c): A few complete solutions which are drawn with the parameters at Ac = 
1.7, an = 0.05, rh = 0.2 and C = 5 and only the inner sonic point is varied: (a) a;j„ — 
2.545, (b) Xin = 2.55 and (c) Xin ~ 2.555. The corresponding shock locations are (a) 
Xs = 48.199, (b) Xs = 27.8854 and (c) Xs = 18.6445 respectively. Vertical dashed hnes 
show the shock transitions. 



parameters are C = 1 (dot-dashed), 5 (long-dashed) and 10 (sohd) respectively. In 
this case, the parameter space shrinks drastically when cooling is enhanced. 



7. Sonic Points for Super-Keplerian Flows? 

It is generally accepted that at the sonic points, the flow must be sub-Keplerian 
(C90b). However, if the cooling is very efficient, this requirement may be violated. 
In Fig. 11, we show that when ( is very high, of the order of a few hundreds, the 
sonic point may have angular momentum above the Kcplcrian distribution (dotted 
curve). Particularly important is that this is possible if the sonic point is located 
near the black hole horizon. The implication of this is not obvious. Does it mean 
that very cold flow can be Kcplcrian or super-Keplerian throughout its journey? If 
so, can it spin up the black hole faster that what is presumed so far? This point is 
to be addressed in future. 
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8. Critical Cooling and Sub-division of the Parameter Space 

We have already indicated that cooling and heating have opposite effects in decid- 
ing the solution topologies, but one does not exactly cancel the other effect. When 
the cooling is enhanced for a given viscosity parameter, the possibility of shock 
formation is eventually reduced. This is shown in Fig. 12. Here the critical cool- 
ing parameter ( = (cri is plotted against the specific angular momentum for two 
different viscosity parameters. Solid curve is for an = 0.01 and the dashed curve 
is for an = 0.05. All possible inner sonic points are considered. The region below 
the curve contains topologies which are closed and therefore standing or oscillating 
shocks could be possible while the region above the curve allows solutions with open 
topologies. We note that for smaller A, the critical cooling factor is small which for 
larger A. This is expected since the possibility of shock formation is enhanced with 
A in general. When a is higher (cri is lower. This indicates that general reduction 
of the parameter space due to higher viscosity^ . is not totally compensated for by 
cooling effects. 
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Fig. lO(a-b): Modification of the parameter space for shock formation when cooling is 

varied, (a) an = 0.01 is chosen and ^ = 0.01, (dot-dashed) 0.1 (long-dashed) and 1 (solid) 
(left panel), (b) ajj = 0.05 is chosen and C = 1 (dot-dashed), 5 (long-dashed) and 10 
(solid) (right panel). 

9. Concluding Remarks 

In this paper, we studied the dissipative accretion flow in presence of viscous heat- 
ing and bremsstrahhmg cooling processes. Viscosity tends to heat the flow, thereby 
reducing the Mach number. Cooling, on the other hand, increases the Mach num- 
ber. Thus, formation of shock, which involves a supersonic to subsonic transition 
is affected by heating and cooling. We classified the parameter space in terms of 
whether three sonic points or shocks can form or not. We discovered a completely 
new topology in which matter, coming from a large distance, is connected to the 
black hole horizon as a normal solution, but it has multiple valued Mach number 
solution. We find that cooling can 'undo' the effect of heating on topological prop- 
erties, only to certain extent. If the viscosity is high enough, then no matter how 
much cooling is used, the parameter space shrinks. We have also found that for a 
given set of flow parameters, there is always a critical cooling factor which separates 
the parameter space into closed regions, one with a closed and the other with an 
open topology. 
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Fig. 11: Example of parameters which produce transonic solutions with super- Keplerian 

angular momentum. The dot-dashed curve is the Keplerian distribution. Solid curves, from 
the bottom to the top, are for = 400, 500 and 600 respectively. 

It is generally believed that only the sub-Keplerian flows can pass through the 
sonic points. However, we find that when the flow is very strongly cooled, even 
super-Keplerian flows can also pass through the inner sonic point. This is a new 
result and may be significant in evolution of the spin of the accretion black holes. 

This work is partly supported by a project (Grant No. SP/S2/K-15/2001) 
funded by Department of Science and Technology (DST). 
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Fig. 12: The variation of the critical coohng parameter as a function of the specific angular 
momentum at the inner sonic point for two different viscosity parameters. Solid curve is 
for an = 0.01 and the dashed curve is for an = 0.05. The region below the curve contains 
topologies which are closed and therefore standing or oscillating shocks could be possible 
while the region above the curve allows solutions with open topologies. 
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